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Abstract. We investigate the effect of external pressure on magnetic order in 
undoped LnFeAsO (Ln = La, Ce, Pr, La) by using muon-spin relaxation measurements 
and ab-initio calculations. Both magnetic transition temperature T m and Fe magnetic 
moment decrease with external pressure. The effect is observed to be lanthanide 
dependent with the strongest response for Ln = La and the weakest for Ln = Sm. The 
trend is qualitatively in agreement with our DFT calculations. The same calculations 
allow us to assign a value of 0.68(2) hb to the Fe moment, obtained from an accurate 
determination of the muon sites. Our data further show that the magnetic lanthanide 
order transitions do not follow the simple trend of Fe, possibly as a consequence of the 
different /-electron overlap. 
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1. Introduction 

Superconductivity in iron based (IB) pnictides and chalcogenides originates from a 
magnetic parent, as in the cuprates. This remarkable analogy hides the well known fact 
that IB materials are always metallic, pointing to a different role of electron correlations 
in the two cases. Another noteworthy difference is that the phase diagram of the IB 
compounds, from the magnetic to the superconducting phase, can be spanned not only 
by doping, as in the case of cuprates, but also by pressure [HE1IB1II1E1E1IH1IH] and by 
nominally isovalent substitutions [101 HU [12] . 

The LnFeAsO i- x F x family where Ln is a lanthanide, also referred to as 1111, 
is particularly relevant because it still holds the record superconducting critical 
temperature T c = 56K [5J. The magnetic and superconducting properties of 1111 
are however the least sensitive to pressure [TJ [HI El US] among IB materials, and the 
pressure dependence of the magnetic moment in the parent members is not known. This 
motivates the present investigation on the effect of pressure on the magnetic ordering 
temperature and on the magnitude of the zero temperature ordered moment for Ln = 
La, Ce, Pr, Sm. 

These four compounds have already been investigated at ambient pressure by 
neutron scattering, providing the ordering vector of the spin density wave (SDW), [1/2 
1/2 1/2] and [1/2 1/2 0] in the tetragonal setting, for La and for Pr, Ce, respectively 
[14] . The value of the magnetic moment on Fe was also obtained and apparently it 
varies between 0.25 /ig and 0.8 //g, but the influence of the rare earth moment on the 
measurements and some earlier results may be questioned. Ambient pressure muon spin 
rotation /zSR |T31 EH EH US] and Moessbauer measurements [211] indicate instead that 
the magnetic moment on iron has very similar values in all compounds (within 0.06% 
of each other). At low temperature (T < 5 K) muons detect also the ordering of the 
rare earths. In the special case of Ce the polarization of the rare earth by the ordered 
Fe lattice, indicating a strong Fe-Ce hybridization, provides an additional contribution 
to the local field at the muon site, observed already well above the magnetic ordering 
of the Ce sublattice. [H] 

In this work we describe the muon results under pressure, from which we extract 
the pressure and the temperature dependence of the Fe magnetic moment, on the basis 
of an ab-initio determination of the muon site. The same ab-initio calculations provide 
insight in the discussion of the Fe moment behaviour. 

2. Methods 

In pressure experiments a large fraction of the muons, roughly 50%, stops in the pressure 
cell (described elsewhere in details [9] [21]), providing a background contribution that 
has to be separated from the sample signal in the data analysis. 

In these conditions the most straightforward method to determine the magnetic 
transition temperature under pressure is by means of /iSR measurements in a weak 
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Figure 1. (Colour on-line) Typical time evolution of the WTF-/1SR asymmetry for 
the LaFeAsO sample with an external field of Bo = 5 mT, applied perpendicular to the 
initial muon polarisation. Above the magnetic transition both the pressure cell and 
the sample contribute to the initial oscillating amplitude, whereas only the pressure 
cell contributes at low temperatures. 



external magnetic field H, applied orthogonal to the initial muon spin direction, in the 
so called weak transverse field (WTF) configuration. In this case the contribution to 
the asymmetry from muons experiencing a vanishing internal spontaneous magnetisation 
can be accurately determined. Muons stopping in a non magnetic environment produce 
long lived oscillations, which reflects the coherent muon precession around the external 
field -B = HqH, whereas muons stopping in magnetically ordered parts of the sample 
give rise to a more complex, distinguishable signal, reflecting the vector combination 
B^ — Bi + Bq of the internal and applied fields. The random orientation of the grains 
in a polycrystalline sample leads to a broad distribution of precession frequencies. 

A distinct experiment can be performed in zero applied field (ZF). In this case 
muons precess around the internal field alone, which may be then more accurately 
extracted from the oscillating muon asymmetry, at the price of a less straightforward 
identification of the transition temperature, where the precessing component disappears. 
Ambient pressure ZF /iSR experiments can be further calibrated in separate experiments 
in the absence of the pressure cell, stopping therefore nearly all the muons in the sample. 

In the following we describe in more details the experimental results for a couple 
of our LnFeAsO samples, namely Ln=La, Ce, and we conclude with a summary of the 
results for the entire series of samples 

2.1. Determination of magnetic transition T m 

Let us first describe the WTF procedure in the representative case of the Ln=La sample. 
We performed WTF /zSR measurements between 5 K and 200 K in an external field of 
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Figure 2. (Colour on-line) Temperature dependence of the magnetic volume fraction 
at different external pressure for the LnFeAsO samples with Ln=La, Ce, Pr, Sm. 



/j,qH = 5 mT, perpendicular to the initial muon spin direction S^. The three panels of 
Fig. [T] display the time evolution of the WTF-/iSR asymmetry at different temperatures, 
for three fixed applied pressure values. The solid curves in the figure represent the global 
fit at each pressure to the following function: 

AvtfW = (a PC e- CT ^ 2 /2 + apM e -wj C os( 7j B t) 

+ a T e- CT ^ 2/2 cos^Bit) + a L e" Ai * (1) 

where , y/2ir = 135.5 MHz/T is the muon gyromagnetic ratio. The first term is the 
contribution of muons stopping in the pressure cell (cell asymmetry), and its small 
relaxation rate, o>c = 0.5 /is -1 , has been calibrated separately with the empty cell. The 
initial asymmetry ape is fitted globally, imposing a common value for all temperatures. 
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Figure 3. (Colour on-line) Pressure dependence of the magnetic transition 
temperatures T m and transition widths AT m as a function of external pressure for 
the LnFeAsO samples with Ln=La, Ce, Pr, Sm 

The second term is due to muons implanted in the part of the sample which is still in the 
paramagnetic phase (paramagnetic asymmetry Opm)- The global procedure converges 
rapidly thanks to the fact that the paramagnetic asymmetry vanishes for T <C T m , as 
it was verified for all the samples at ambient pressure in a separate WTF experiment 
without the pressure cell. Muons stopping in a local environment that is magnetically 
ordered are described by the last two terms (the magnetic asymmetry), which reflect 
respectively the transverse (B^ JL S^) and the longitudinal (B^ \\ S^) components of 
the total local field, averaged over all grain orientations. 

The temperature dependent magnetic fraction of the sample is evaluated as 
1 — clpm/uo where a is the total muon asymmetry from the sample, obtained from 
the global fit at high temperatures, well above the magnetic transition T m . The same 
procedure has been applied to all our samples. The values of the magnetic volume 
fraction for LnFeAsO Ln=La, Ce, Pr, Sm are displayed as a function of temperature 
in Fig. [2] for the different values of the applied pressure. The magnetic transitions and 
their widths have been evaluated by the best fit of each experimental data set to the 
function erf[(T — T m ) / (\/2AT m ))], which accounts for a Gaussian distribution of T m . 
These parameters are displayed in Fig. |3j together with their linear regressions (solid 
lines for T m ). The dashed line for AT m is the pressure dependence of the width of the 
transition and it is indicative perhaps of the increasing inhomogeneity of the pressure 
within the cell. 
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Figure 4. (Colour on-line) Time evolution of the ZF-/1SR asymmetry for the LaFeAsO 
sample at ambient (left) and maximum (right) pressure, at three representative 
temperatures. 



2.2. Determination of the internal field at the muon sites Bi 

The most accurate determination of the internal fields and how they change as a 
function of temperature and pressure is obtained by separate ZF /xSR experiments. 
The spontaneous internal fields Bi at the muon site are proportional to the intensity 
of the staggered magnetic moment on iron, m, with an additional contribution from 
magnetic rare earths below their ordering temperatures. Representative ZF-/iSR spectra 
are shown in Fig. [4] and Fig. |5j for Ln=La and Ce respectively. Equation [T] may be 
rewritten now for B = as: 



A ZF (t) = a PC e- a rc t2 / 2 + a PM e" w 

+ a " e- CT «* 2/2 cos( 7 5 4 t) + a L e~ XLt (2) 

i=l,2 

where higher statistics ZF experiment allows us to recognise two distinct transverse 
components in the magnetic asymmetry, indicated by subscripts i — 1, 2. They are due 
to the two known muon stopping sites in 1111 compounds [IB] , with B\ corresponding 
to a site within the FeAs layer and B2 corresponding to a site close to O ions, for which 
we find an occupancy ratio of aTi/a T 2 ~ 4. The corresponding longitudinal components 
cannot be resolved in the last term of the Eq. |2j 

The low field B 2 is accurately detected only for the Ln=La sample and for T < 70 
K (Fig. El), reaching the value -62(0) ~ 18 111T at ambient pressure, in agreement with 
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Figure 5. (Colour on-line) Time evolution of the ZF-^SR asymmetry for the CeFeAsO 
sample at ambient (left) and maximum (right) pressure, at three representative 
temperatures. 



Ref. p2]. In the other investigated LnFeAsO compounds this low precession is over- 
damped due to the broader field distribution at this muon site. This is justified by the 
stronger magnetic perturbation provided by the magnetic rare earth, much closer to this 
than to the other site. For the same reason the muon field B\ can be easily detected 
almost up to the transition temperature. The two fields are displayed as a function of 
temperature in Fig. [6j 

We are interested in particular in the effect of pressure on the zero temperature 
value of the staggered magnetic moment of Fe, m oc Bi(0). The extrapolation to zero 
temperature must avoid including any very low temperature variation of the internal 
field, due to the onset of rare-earth ordering. To this end the internal field has been 
fitted for T > Tjf 1 to the phenomenological fit curve 



B i (T)=B i (0) [1-{T/T, 



al/3 
mj J i 



(3) 



that interpolates between a low temperature Bloch-like regime and the critical regime, 
close to T m . 

In the case Ln = Ce the molecular field of Fe on the paramagnetic Ce moments [TS] 
produces and additional non negligible Curie- Weiss contribution from Ce at the muon: 



Bi{T) = Bi(0) 



1 - 
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(4) 
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Figure 6. (Colour on-line) Temperature dependence of the two internal fields, Bi, 
i = 1,2: open symbols from ZF fiSK, filled symbols from WTF /iSR. Solid lines are 
the phenomenological fit of Bi(T) to Eq. [3] (Ln = La, Pr, Sm) and to [4] (Ln = Ce); 
the dashed line for Ln = La is the fit of the lower B 2 (T) field. 



The fits have been done by using the same T m values reported in Figure [3] The power 
coefficients are 1.5 < a < 2.5 and 0.15 < (3 < 0.2 for all the samples. The results of 
these fits are shown in Figure [7] for -Bi(O) as a function of pressure for the all samples 
under investigation. 

2.3. Determination of the magnetic transition of Ln=Sm,Pr,Ce sub-lattice, Tjf 1 

It is known that the at ambient pressure the magnetic Ln ions orders antiferromagnet- 
ically below Tn = 4.6, 4.4 and 11 K for Ln=Sm, Pr and Ce respectively [15] . This 
magnetic transition introduces a Ln-dependent contribution to the muon precession fre- 
quency in the ZF-/zSR signal, i.e. to the internal field of Eq. [2j At ambient pressure 
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Figure 7. (Colour on-line) Pressure dependence of the internal field at the most 
populated muon site, extrapolated at zero temperature, -Bi(O), for the LnFeAsO 
samples with Ln=La, Ce, Pr, Sm. Straight lines are the best fit to a linear behaviour. 



the magnetic transitions of Ln = Pr, Ce are second-order-like [18] and produces a re- 
duction of B^T) in the case of Pr, whereas for Ce the field shows a broad peak. Panel 
a and b of Figure [8] show that the same behaviour is found also under external pressure. 
T^ n is found to increase for both rare earths, being T^ r = 12.1(1) and T^ e = 6.2(2) K 
at the maximum applied pressure p = 2.3 and 2 GPa, respectively. 

On the other hand, the magnetic transition of Sm in SmFeAsO at ambient pressure 
is discontinuous at T^ m (first-order) [TH] and it is accompanied by a change of the 
magnetic unit cell. Below T^ m the single muon precession frequency line due to the 
order of Fe moments develops two smaller amplitude satellite lines. Unfortunately 
these satellites cannot be resolved in the presence of the large pressure cell background 
signal, and show up just as an increased static relaxation. The sharp change of the 
depolarisation rate of the transverse component, proportional to the width of the of 
the internal field distribution, is however sufficient to signal the onset of Sm order and 
its pressure dependence. Panel a in Fig. [8] shows that the maximum applied pressure 
p = 2.3 GPa produces in this case a reduction of the magnetic transition to T^ m = 3.8(1) 
K. 

3. Electronic structure and muon field calculations 

Electronic structure calculations are performed within the Density Functional Theory 
(DFT) and local spin density approximation as parametrised by Perdew and Wang [22] 
plus Hubbard corrections (LSDA+U) [23] as implemented in the all-electron LAPW 
code Wien2K [21]. The muffin-tin radii for Ln = La, Ce, Sm, for Fe, As, and O are 



Pressure effects in LnFeAsO 



10 



2600 

2200 

1800 
1700 
1500 
1300 
1100 
50 



-A-O.OGPa 
° 2.0 GPa 



CeFeAsO - 




(a) 

i i i i i i i I i i i I i i 1 1 1 1 i 



-0.0 GPa 
-1.1 GPa 
-2.3 GPa 



I 

PrFeAsO ■ 



V r =11.6(1)K 

-a—a—°kx, o- 




T(K) 



Figure 8. (Colour on-line) Temperature dependence of the internal He\d,B^(T), for 
Ln = Ce (a) and Pr (b) and of the depolarisation rate of the transverse fraction, ar(T), 
for Ln = Sm (c) nearby the antiferromagnetic transition of the Ln sub-lattice Tjf 1 



chosen equal to 2.3, 2.2, 2.0, and 1.9 Bohr respectively. -Rmt x kmax = 7 is used as the 
plane-wave cutoff. The U on-site Coulomb parameter is set to 6 eV for Ce and Sm, and 
the collinear magnetic formalism is adopted. 

The DFT calculation allows a straightforward identification of the muon site, 
corresponding in first approximation to the minimum of the electrostatic potential (i.e. 
the sum of the Hartree potential plus the local pseudopotential), reversed in sign, since 
the muon has a positive charge. In the case of the 1111 structure we identify three 
minima, the lowest two replicating those found in Ref. [TS] with a simple Thomas-Fermi 
approach (we quote here results for Ln=La). 

Given the small muon mass (m = 105.7 MeV/c 2 ), we investigate possible effect of 
the zero-point-motion on the stability of the identified minima. We model the potential 
around each minimum as an anisotropic harmonic well V(r) = \m{oj 2 x 2 +u 2 y 2 +uj 2 z 2 ) + 
Vq with eigenvalues given by E(n x , n y , n z ) = h[u x {n x + 1/2) + to y (n y + 1/2) + u z {n z + 
1/2)] + Vq. The location, the minimum potential Vo, the zero point energy Eq and 
the first excitation energy E\ are shown in Tab. [TJ The actual values of E\ indicate 
that in every site only the ground state is occupied at room temperature. From these 
numbers we estimate the extent of the muon wave-function around each site by drawing 
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Figure 9. (Colour on-line) Muons sites in the LaFeAsO unit cell by DFT. The golden 
light shaded areas represent the volume of zero-point displacement for the muon around 
the minimum. 

the isosurface corresponding to V(r) = V + E . Figure [9] shows these isosurfaces as a 
light-shaded areas. Assuming that in site a the muon overtakes any barrier lower or 
equal to Eq , we see that site A and B are disconnected, while site C and A form an 
interconnected network where the migration barrier is lower than Eq' . Therefore site 
C is unstable and we are lead to conclude that, while the muon potential has three 
nonequivalent minima, only site A and B are surrounded by barriers high enough to 
guarantee muon confinement. These sites correspond roughly to those identified in 
Ref. [IE]. 

Table 1. Muon sites: cell coordinates x, y, z, potential V (eV) referred to the 
minimum, zero point energy, first excited level Eq,E\ (eV), and local dipolar field, 
assuming a value of the Fe moment fi — 0.68 (is (see text). 
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0.75 
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0.63 
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0.25 


0.25 
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0.55 


0.51 


0.09 


29 


He 


0.5 


0.0 


0.703 


0.94 


0.54 


0.25 





If we now calculate the internal fields, Bi as the distant-dipole contribution from 
the Fe moments in the known [13] magnetic structures, a Fe staggered magnetic 
moment m = 0.68(2) fis along (110)r is required to obtain the experimental value 
-Bi(O) = 0.165(5) T at site A. By the same token site B cannot be the main muon site, 
since it requires an unacceptably high Fe moment of 3.6 yU_e to reproduce -Bi(O). The 
field reported in Tab. [T]for site B, with 0.68 /is on Fe, agrees within a factor 1.6 with 
B 2 (0). 
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Figure 10. (Colour on-line) LaFeAsO band structure under pressure. Continuous 
(red) and dashed (black) lines show the band structure for the 2 GPa compressed 
structure and equilibrium one respectively (Colour on-line). 

We are further interested in the variation of the magnetic properties with pressure. 
Since the experimental values for both lattice parameters and internal position are not 
available for all of the compounds under investigation, theoretical position are used 
instead, also for the ambient pressure condition. The Wien2K code allows for the 
optimisation of the internal degrees of freedom by the use of the Hellmann-Feymann 
forces. For each compound calculation of the total energy and optimised internal 
structure was performed for different values of the lattice parameters a and c for the 
tetragonal structure (we neglected the orthorhombic distortion). Upon interpolation of 
the total energies we obtain the optimised structures at different volumes. The volume 
versus pressure dependence is obtained by a fit of the Murnaghan equation of state, 
which yields the theoretical structural parameters. 

We find that pressure reduces the magnetic moment on Fe in all three materials 
under investigation, confirming previous calculations for RE=Ln [25 1 126 | 127]. This is in 
agreement with the Stoner criterion: the energy gain for the magnetic state correlates 
inversely with the width of the d-bands, up to a point where the system reverts to the 
paramagnetic state. Indeed, in pure Fe the reduction of volume induced by applied 
pressure increases the width of the rf-bands. In the present case, the maximum 2.3 
GPa pressure is not sufficient to induce a paramagnetic transition, but the effect is in 
agreement with this qualitative argument. 

This picture is confirmed in the case of LaFeAsO by the comparison of the band 
structures for the equilibrium structure (solid) and under 2 GPa pressure (dashed), 
shown in Fig. [3] We see that the width of the valence manifold increases about the F 
and A point as evidenced by the vertical arrows. 
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Figure 11. (Colour on-line) Pressure dependence of magnetic temperature of the 
FeAs layers, dT m /dp, and staggered magnetisation at zero temperature, dM s /dp, rates. 
Lines are guides to the eye. 

The ab-initio muon site identification agrees with previous estimates j 18 . yielding 
the same determination of the iron moment, compatible with the upper end of the 
neutron diffraction values [14J. Since muons conclusively determine that all Ln 1111 
compounds have very close internal fields our calculations guarantee that they must 
have correspondingly close values of the iron local moment. 

Our determination of the iron moment is obtained by comparing the calculated 
distant dipole field with the experimental value from /zSR. In principle two additional 
contributions are present: the contact Fermi field, written as Bp = 2fioTrip e (p^(0) — 
Pl(0))/3 for collinear magnetic structures, where p^i( r ) is the up(down) spin density at 
a position r and the muon site is at r = 0; the anisotropic on-site pseudo-dipolar field, 
i.e. the the dipolar field averaged over the central muffin tin. We notice, however, that 
the three sites listed in Tab. [T] share high x, y symmetry. Thanks to this condition the 
experimental staggered collinear structures yield Bp = by symmetry in all of them, 
fully justifying our neglect of this contribution. Furthermore the DFT calculations 
shows that the p-wave contribution at the muon, dominant in the pseudo-dipolar field, 
is also negligible. Consequently the distant dipole field is a very good estimate of Bi and 
within this approximation /zSR yields a direct measurement of the moment on iron. We 
note that its value, mp e = 0.68(2) ps, is in very close agreement with the most recent 
neutron scattering [28J, mp e = 0.63(1) pb, and NMR [29J, m Fe = 0.58(9) ps- 

A more precise calculation, beyond the scope of the present work, involves the 
expectation value of the total local field over the muon ground state, roughly represented 
in Fig. [9] by the shaded volumes around each site. This is a second order effect that 
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is certainly more significant at site B, in view of it larger and non isotropic zero point 
volume, compared to that of site A. As we pointed out the dipolar field at B site (Tab. [T]) 
overestimates -82(0) by a factor 1.6. However, the average field at the boundary of the 
zero point volume in B is reduced by a factor 0.5, hence the accurate inclusion of the 
muon zero point motion will certainly improve the agreement between calculations and 
experiments also for this site. 

The results of our study show that for all the compounds under investigation, 
LnFeAsO with Ln = La, Ce, Pr and Sm, both the staggered magnetisation, M s , and the 
magnetic transition temperature of the FeAs layers, T m , are progressively diminished as 
a function of the external pressure. Up to the maximum applied pressure of p m 2.3 GPa 
all the samples display only a modest variation of the magnetic properties, indicating 
that the complete destruction of magnetism requires many GPa. By DFT calculations 
we have shown that the reduction of the Fe magnetic moment as a function of pressure 
is mainly related to the broadening of the ci-bands. As expected in an itinerant magnet 
such as LnFeAsO, this increases the kinetic energy and reduces the exchange energy for 
the spin alignment. 

We experimentally found that the rate of the reduction of both the magnetic 
temperature, dT m /dp, and the staggered magnetisation, dM s /dp, are Ln dependent. 
Figure 11 shows that both increase almost linearly as a function of the ionic radius, 
i.e. of the chemical pressure which increases from La to Sm. This behaviour indicates 
that the chemical pressure makes the structure of the latter more rigid, reducing the 
effect of the external pressure. A similar effect has been recently reported for Ce alone 
[30J, where the magnetic transition is detected indirectly through the peak in resistivity 
at T*, associated with the structural transition, precursor of T^. The shift in T* from 
145(5) to 124(5) K under the application of 2.2 GPa agrees well with that reported by us 
in Fig. [3j from 127 to 110 under 2.0 GPa. It is interesting to note that the the sensitivity 
of magnetism to pressure correlates inversely with the optimal superconducting T c of 
each Ln family. 

Our results show that the effect of the external pressure on the Neel transition 
temperature of the magnetic Ln sub-lattice is also Ln dependent displaying an increase 
of Tjj n for both Ce and Pr and a decrease for Sm. Figure [8] shows that the biggest 
variation is observed for Ln=Ce, with a rate of dT^ e /dp = 0.8(1) K/GPa, in agreement 
again with the transport measurements reported in Ref. [30]. This variation has been 
related to the increase of the f-d hybridization between the Fe and Ce orbitals, which 
leads to the enhancement of the indirect exchange coupling among the Ce ions [30]. The 
same effect in the Pr and Sm cases is expected to be significantly smaller since the degree 
of f-d hybridisation is much lower [31]. This can explain the small increase observed 
for Ln=Pr in Fig. [8Jd, which leads to dT^ r /dp = 0.22(5) K/GPa. On the other hand 
Ln = Sm seems to undergo a first order magnetic transition, according to Ref. [18] . If 
this is the case it should be accompanied by a spontaneous deformation of the crystal 
structure with a concomitant variation of the volume cell. The external pressure would 
then contrast the increase of the volume cell and hence lower the associated T^ m , which 
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would explain the observed negative dT§ m /dp = —0.3(1) K/GPa. The occurrence of a 
change of the volume cell deserves an experimental confirmation. 

In conclusion, we have measured the evolution of the magnetic properties of both 
the FeAs and LnO layers and performed DFT calculations for the Fe magnetic moment 
and the muon sites in LnFeAsO (Ln = La, Ce, Pr, Sm) as a function of external pressure 
up to ~ 2.3 GPa. Our results show that the behavior of all the magnetic properties 
under pressure is Ln dependent. Regarding the FeAs layers, both the magnetic transition 
temperature T m and the Fe magnetic moment decrease as a function of pressure for all 
the Ln, but with a pressure rate which scales almost linearly with the Ln ionic size, being 
maximum for the larger Ln = La atoms. On the other hand, the pressure dependence 
of the antiferromagnetic transition temperature of the LnO sublattice Tjj n for magnetic 
Ln = Ce, Pr, Sm ions, do not follow the simple trend of Fe, but it seems to be sensitive 
to the different f-electron overlap. In addition Ln = Sm displays a neg ative dT§ m /dp 
rate which might be inferred to the first-order nature of the sublattice Sm transition. 
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